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ABSTRACT 

Context. 

Aims. We present a study of the molecular gas and interstellar dust distribution in the environs of the Hii region NGC 3503 associated with the open 
cluster Pis 17 with the aim of investigating the spatial distribution of the molecular gas linked to the nebula and achieving a better understanding 
of the interaction of the nebula and Pis 17 with their molecular environment. 

Methods. We based our study in '^CO(l-O) observations of a region of ~ 0?6 in size obtained with the 4-m NANTEN telescope, unpublished radio 
continuum data at 4800 and 8640 MHz obtained with the ATCA telescope, radio continuum data at 843 MHz obtained from SUMSS, and available 
IRAS, MSX, IRAC-GLIMPSE, and MIPSGAL images. 

Results. We found a molecular cloud (Component 1) having a mean velocity of -24.7 km s"' , compatible with the velocity of the ionized 
gas, which is associated with the nebula and its surroundings. Adopting a distance of 2.9 ± 0.4 kpc the total molecular mass and density yield 
(7.6 + 2.1) X 10^ Mo and 400 ± 240 cm"\ respectively. 

The radio continuum data confirm the existence of an electron density gradient in NGC 3503. The IR emission shows the presence of a PDR 
bordering the higher density regions of the nebula. The spatial distribution of the CO emission shows that the nebula coincides with a molecular 
clump, with the strongest CO emission peak located close to the higher electron density region. The more negative velocities of the molecular gas 
(about -27 km s"' ), is coincident with NGC 3503. Candidate YSOs were detected towards the Hii region, suggesting that embedded star formation 
may be occurring in the neighbourhood of the nebula. The presence of a clear electron density gradient, along with the spatial distribution of the 
molecular gas and PAHs in the region indicates that NGC 3503 is a blister-type Hii region that probably has undergone a champagne phase. 
Conclusions. 

Key words. ISM: molecules. Radio Continuum: ISM, Infrared: ISM, ISM: Hn regions, ISM:individual object: NGC 3503 , Stars: Pismis 17 
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It is well established that massive stars have a significant impact 
, on the dynamics and energetics of the interstellar medium (ISM) 

■ surrounding them. In the classical scenario, OB stars are born 
' deeply buried within dense molecular clumps and emit a copi- 
ous amount of far ultraviolet (FUV) radiation (hv > 13.6 eV). 
FU V photons ionize the neutral hydrogen, creating an Hii region 
which expands into the molecular cloud due to the pressure dif- 
ference between the molecular and ionized gas. At the interface 
between the ionized and molecular gas, where photons of lower 

[ energies (6 eV < hv < 13.6 eV) do minate, photodissociation re- 
. gions (PDR) are known to exist (see lHollenbach & Tielensiri997l 
for a complete review). 

The molecular compo nent of Hii region s has been studied 

■ in d etail by many authors ( Keto & Holll989l : iRusseil & CastetsI 
i 12004 ). Many efforts have been devoted to observing the interac- 
tion between the Hii regions and their parental molecular clouds 
and some evidence has been presented showing that molecular 
gas near Hn regions may be kinematically disturbed by several 
km s~^ (Elmegreen & Wand ll987l) . Understanding the complex 
interaction between massive young OB stars, Hii regions, and 
molecular gas is crucial to the study of massive star formation 
and the impact of massive stars on their environment. 

NGC 3503 (=Hf 44=BBW 335) (l,b = 289?51, +0^12) is a 
bright and small (~ 3' in diameter) optical emission nebula 
(iDrever & Sinnottlll988 ^ ionized by early B-type stars belong- 
ing t o the open cluster Pis 17 (Herbst 1975; Pinheiro et aL, 
12010 ) and it is seen projected onto an extended region of 
medium brightness in Ha known as RCW54. This last re- 



gion is centered at (I, b) - (289°4, -0^6), and has dimensions 
(Al X Ab) = (3°5 X r.O) (Rodgers et al. 1960). The point source 
IRAS 10591-5934, which appears projected onto NGC 3503 , is 
the infrared counterpart of the optical nebula. For the sake of 
clarity, we show the main components and their relative location 
in Fig. [1] superimposed on a DSSR image of the region of inter- 
est. The dashed white square delimits the region studied in this 
paper 

In their Southern Hemisphere catalogue of bright- 
rimmed clouds (BRCs) associated with IRAS point sources, 
Sugitani & Ogura ( 1994) quote SFO 62 as related to NGC 3503. 
BRCs are defined as isolated molecular clouds located at 
the edges of evolved Hii regions, and are suspected to be 
potential sites of star form ation through the r adiation-driven 
implosion (RDI) process ( Sugitani e t al.l Il99lt 'Urauhart et alj 
2009). Making use o f the Australian Telescope Compact Array 
(ATCA), iThompson etafl im)40 caiTied out radio continuum 
observ ations towards all BRCs catalogued by S ugitani & Ogural 
11991 . These authors classified SFO 62 as a broken-rimmed 
cloud associated with an evolved stellar cluster that is about 
to disr upt its natal molecular cloud. Latter on, Urauhart et aT| 
dloosi) claimed that SF062 was incoiTectly classified as a BRC, 
and that the RDI process is not working in this region. 

Using the [Snl i6716/A6731 Hne ratio of NGC 3503, 
ICopetti et al.l (l2000h detected a significant electron density de- 
pendence on position, which according to the authors, could 
be interpreted as a radial gradient. The authors claimed that 
NGC 3503 may be a candidate of showing a "champagne 
flow". In their Ha survey of the Milky Way towards / - 290°, 
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Fig. 1. DSSR image of the brightest part of RCW 54. The image is ~ 1?5 x l^S in size, centered at {l,b) = (289?5, 0?25). Our region 
of interest is delimited by the white dashed square. 



iGeorgehn et al.l (l2000l) reported a radial velocit}Q of -2 1 km s ' 
towards NGC 3503 , and proposed that NGC 3503 is hnked to a 
complex of Hii regions placed at a distance of about 2.7 kpc, 
with radial velocities of about -25 km s ' . Published distance 
dete rminations of NGC 3503 vary between a minimum of 2.6 
kpc (|Herbsllll975h and a maximum of 4.2 kpc dMoffat & Vogtl 

: 0.4 kpc 



9751). In this paper we shall adopt a distance of 2.9 : 
jPinheiroet al.i2010l) . 



NANTEN ^^CO (J=1^0) (HPBW ^ 2.7 line observations 
were carried out by lYamapuchi et al.l (Il999h towards 23 Hn re- 
gions associated with 43 BRCs catalogued by Sugitani & Ogura| 
(|1994|) . with a grid spacing of ~ 4' x 4'. The aim of that work 
was to investigate statistically the dynamical effects of Hii re- 
gions on their associated molecular clouds and star formation. 
A square region of ~ 0°4 in size towards NGC 3503 was ob- 
served and a single and slightly extended 6' diameter molec- 
ular cloud having a peak radial velocity of -24.9 km s"' was 
found. Adopting a distance of 2.9 kpc, this molecular structure 
has a linear radious of 2 pc and a total mass of ~ 500 . 
The authors rejected this cloud from the analysis sample, since 



Radial velocities are referred to the Local Standard of Rest (LSR) 



it was detected at only one position. '^CO, '^CO, and '^CO 
(J=l— > 0) MOPRA po s ition -switched observations were carried 
out bv lUrquhart et al.l (12009) . The ON position was centred on 
the position of the IRAS source associated with NGC 3503. The 
least abundant of the isotopes was not detected, whilst the emis- 
sion lines of the other two isotopes show a double-peaked pro- 
file. From a Gaussian fitting, their peak radial velocities are ~ 
+20kms"' and -25.6 km s"' . 

Molecular line observations provide an invaluable sup- 
port in pursuing a better understanding of the interaction of 
Hn regions with their s urrou ndings. Previou s observations of 
lYamaguchi et al.1 d 19991) and lUrquhart etal] [2009.) . although 
providing important information about the molecular gas asso- 
ciated with NGC 3503 , do not offer a complete picture of the 
molecular environment of the nebula. Bearing this in mind, the 
goals of this paper are twofold, namely: aj to map the spatial dis- 
tribution of the molecular gas associated with NGC 3503 and to 
study its physical characteristics, and b) to achieve a better un- 
derstanding of the interaction of NGC 3503 and Pis 17 with their 
molecular environment. These aims were addressed by analyz- 
ing new '^CO (J=1^0) data gathered by using the NANTEN 
telescope to observe a square region ~ 0°7 in size, centred on 
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(l,b) - (289°36, +(f.02). Molecular observations were combined 
with unpublished radio continuum data, and optical, mid, and 
far-infrared archival data. 

The structure of the paper is the following. In Sect. 2, the CO 
and radio continuum observations are briefly described along 
with a short mention to the archival data used in this work. 
The main observational results of the molecular, ionized and 
dust emission are detailed in Sect. 3, while both the comparison 
among these interstellar phases and the action of the stellar ion- 
izing sources of NGC 3503 over its interstellar environs, as well 
as the presence of star formation activity are exposed in Sect. 4. 
Finally, a possible scenario is put forward to explain the stellar 
and interstellar interactions. A summaiy is presented in Sect. 5. 

2. Observations and data reductions 

The databases used in this work are: 

1. Intermediate angular resolution, medium sensitivity, and 
high-velocity resolution '^CO (J=l— »0) data obtained 
wifli the 4-m NANTEN millimeter-wave telescope of 
Nagoya University. At the time the authors carried out 
the observations, April 2001, this telescope was installed 
at Las Campanas Observatory, Chile. The half-power 
beamwidth and the system temperature, including the 
atmospheric contribution towards the zenith, were 2. '7 
(~ 2.3 pc at 2.9 kpc) and ~ 220 K (SSB) at 115 GHz, 
respectively. The data were gathered using the position 
switching mode. Observations of points devoid of CO emis- 
sion were interspersed among the program positions. The 
coordinates of these points were retrieved from a database 
that was kindly made available to us by the NANTEN 
staff. The spectrometer used was an acusto-optical with 
2048 channels providing a velocity resolution of ~ 0.055 
kms"^ For intensity calibrations, a room-temperature 
chopper wheel was employed (Penzias & Burrus 1973|. 
An absolute intensity calibration (Ulich & HaasI ll976t 
Ikutner & Ulich 1981) was achieved by observing Orion KL 
(RA(1950.0) = 5'' 32'"47^0,Dec (1950.0) = -5° 24' 21"), 
and p Oph East (RA(1950.0) = 16'' 29"' 2ff.9, 
Dec (1950.0) = -24° 22' 13"). The absolute radiation 
temperature, T^, of Orion KL and p Oph East, as observed 
by the NANTEN radiotel escope were assumed to be 65 K 
and 15 K, respectively jMoriguchi et al.l l200lh . The CO 
observations covered a region (Al x Ab) of 35 '1 x 35.' 1 
centred at (/, b) = (289^36, +0°.Q2 ) and the observed grid 
consists of points located every one beam apart. A total of 
169 positions were observed. Typically, the integration time 
per point was 16s resulting in an rms noise of ~ 0.3 K. A 
second order degree polynomial was substracted from the 
observations to account for instrumental baseline effects. 
The spectra were reduced using CLASS software (GILDAS 
working group). 

2. Unpublished radio continuum data at 4800 and 8640 
MHz which were kindly provided by J. S. Urquhart and 
M. A. Thompson. The images were obtained in March 2005 
with the Australia Telescope Compact Array (ATCA) with 
synthesized beams and rms noises of 23'.'55 x 18'.'62 and 
0.82mJybeam-i at 4800 MHz and 14'.'73 x 1 1'.'74 and 
0.56 mJy beam"' at 8640 MHz. Radio continuum archival 
images were also obtained from the Sydney Un iversity 
Molonglo Sky Survey (SUMSSfl (iBock et al.lll999h at 843 
MHz, with angular resolution of 45"x45"cosec(5). 

^ ,http://www. astrop.phy sics.myd.edu.au/cgi-bin/postage.pl\ 



3. Infrared data retrieve d from t h e Mid course Space 
Experiment (MSX)^ (Price etdj 1200 ll) . high reso- 
lution IRAS imag es (HIRES D at 60 and 100 //m 
dFowler & Aumannl 1199 4). Spitzer images at 8.0 and 
4.5 jjm from the Galactic Lega cy Infrared Mid-Plane S urvey 
Extraordinaire (GLIMPSE^Q jBeniamin et alj l2003h . and 
Multiband Imaging Photometer for S pitzer (MIPS) images 
at 24 and 70 i m\ from the MIPS I nner Galactic Plane Survey 
(MIPSGALfl (iCarev etal.ll200l . 

4. Optical data r etrieved from the 2nd Digitized Sky Survey 
(red plated ( McLean et al.ll2000l) . 

3. Results and analysis 

3.1. Carbon Monoxide 

3.1 .1 . Individual velocity components 

In order to illustrate in broad terms the molecular structures de- 
tected towards the region under study, a series of CO profiles are 
displayed in Fig.|2] 

Profiles a) and c) show the CO emission along the line of 
sight to NGC 3503 and the bright edge of SFO 62 , respectively. 
In both spectra the bulk of the molecular emission is detected 
between -30 and -20 km s"' , and -Hi 5 to -h25 km s"' , in good 
agreement with previous results bv lUrquharteTan (^009). 

The medium brightness optical region northeast of 
NGC 3503 (from here onwards MBO) though also depicts a 
double peak structure (see profile b), displays some small scale 
structure in the feature peaking at ~ -1-20 km s"' , while the most 
negative CO feature is detected at — 18 km s"' . 

The molecular gas along regions of high optical absorption is 
shown in profiles d) and e). The former resembles the CO spec- 
trum observed towards NGC 3503 , and is characteristic of the 
molecular emission in the region of high optical absorption seen 
northwest of NGC 3503 . The roundish patch of high absorption 
seen at (l,b) = (289?55,-h0°07) only shows the CO peak at posi- 
tive velocities. 

The spatial distribution of the molecular gas observed in the 
three velocity intervals mentioned above is shown in the left 
side panels of Fig. [3] In order of increasing radial velocity, the 
CO components will be referred to as Component 1 (peaking 
at — 25 km s ' ), Component 2 (peaking at ~ -16 km s ' ), and 
Component 3 (peaking at ~ -1-20 km s"' ). To facilitate the com- 
parison between the spatial distribution of molecular and ionised 
gas, the right panels of Fig. |3]show the same molecular contours 
as the left panels superimposed to the DSSR image in grayscale. 
The difference in spatial distribution among the three compo- 
nents is readily appreciated. The molecular gas in the veloc- 
ity interval -29 to -20 km s"' (Component 1) shows two well 
developed concentrations, whose emission peaks are located at 
f/,fej = (289°47,-h0°12) (clump A) and f/>j = (289?32',+0?03) 
(clump B). Component 1 has a very good morphological cor- 
respondence with the high optical obscuration region seen north 
and northwest of NGC 3503 . 

Component 2 peaks at = (289?55, -h0?17) near both 
NGC 3503 and MBO, and is projected onto a region without 
apreciable optical absorption. Finally, Component 3, peaking 
at (l,b) - (289?45, -0?18), does not show a clear morphological 

^ http://irsa. ipac. caltech. edu/ Missions/msx. html\ 
* http://irsa.ipac.caltech.edu/applications/IRAS/IGA/ 
^ http://sha. ipac. caltech. edu/applications/Spitzer/SHA// 
^ http://sha. ipac. caltech. edu/applications/Spitzer/SHA// 
' ^ttp://skyview. gsfc. nasa. gov/cgi -bin/qu ery. pl\ 
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Fig. 2. CO emission profiles towards five positions around NGC 3503. The CO profiles are averaged over a square area ~ 3' in size, 
centered on the black/white dots drawn on the DSSR image (center). The intensities are given as absolute radiation temperature T^. 



correspondence either with ionized regions or areas displaying 
strong absorption. 

For the three molecular components, a mean radial velocity 
(V) weighted by line temperature, was derived by means of 



V 



Tji Tpeaki X Vpeaki 
Tji Tpeak, 



(1) 



where Tpeaki and Vpeaki are the peak T^temperature and the peak 
radial velocity of the /-spectrum observed within the 3 rms con- 
tour line defining the outer border Mean radial velocities of 
-24.7 km s"' , -16.6 km s"' , and -1-19.5 km s"' were obtained 
for Component 1, Component 2, and Component 3, respectively. 

Based on V, we tried to estimate the kinematical distance 
of Componen t 1. Th e analytical fit to the rotation curve by 
iBrand & Blit3 ( Il993h along / = 289?5 shows that radial veloc- 
ities more negative than -15.5 km s ' are forbidden for this 
galactic longitude. However, the ionised gas along this galac- 
tic longitude exhibits velocity departures more negative than -7 
km s"' relative to the circular rotation model, as the cases of the 
Hii regions Gum 37 (l,b) = (290?65, +(f.26) id = 2.7 kpc) and 
Gum38a (l,b) = (291?28, -hO;71) {d = 2.8 kpc), whose radial ve- 
locities are in the ran ge -25 to -28 km s~^ and -25 to - 29 
km s'^ , respectively dGeorgelin et alJ I2OOOI) . Bearing in mind 
that close to the tangential point along / = 290?5 the radial ve- 
locity gradient is very small and that the Hii regions located at 
distances between 2.7 and 2.8 kpc exhibit forbidden radial veloc- 
ities similar to those found for Com ponent 1, we adopt for the 
latter a distance of 2.9 ± 0.4 kpc (iPinheiro et all 120 IQ) . which 



is also the kinematical distance corresponding to the tangential 
point. 

Component 1 shows an excellent morphological resem- 
blance with a region of high optical absorption next to 
NGC 3503 (see Fig. [3]) and its mean velocity is in good agree- 
ment with the velocity of the Ha line towards the nebula 
(-21 km s ' ; iGeorgelin et al.l I2OOO ). It is also worth not- 
ing that the velocity interval of Component 1 is in ex- 
cellent agreement with the velocity of SF0 62 found by 
Sugitani&02ura(1994)J-28< vlsr < 



20kms-i,at '^CO), 



by Yamaguchi et al. (1999) (-24.9 km s"' at '^CO), and by 
Ur quhart et al. (2009) (-25.6 and -25.7 kms-\ at '^CQ 
and '■'CO, respectively). Based on above, we conclude that 
Component 1 is associated with NGC 3503 and its environs. The 
bright rim of SFO 62 clearly follows the southernmost border of 
clump A, which very likely indicates that the molecular gas in 
the southern border of this clump is being ionized, originating 
the bright optical rim. 

On the other hand, considering that the mean radial velocity 
of Component 2 (-16.5 km s"' ) is close to the radial velocity of 
the tangential point at / = 289^5 (-15.5 km s"' ), we suggest that 
Component 2 may be located in the neighborhood of the tangen- 
tial point, close to Component 1 . From the mean radial velocity 
of Component 3, a kinematic distance of ~ 8 kpc is determined, 
indicating that Component 3 is unrelated to NGC 3503 . Very 
likely. Component 3 is associated with the co mplex of Hii re- 
gions at a velocity of ~ +20 km s ' reported bv lGeorgelin et al.) 
(I2000h . 
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Fig. 3. Upper panels (Left): Averaged T^in the velocity range 29 to - 20 km s"' (Component 1). Contour levels start at 0.42 

K (23 rms) and the contour spacing is 0.6 K. The beam size of the CO observations is shown by a circle in the lower left corner. 
(Right): Overlay of the revalues in the same velocity interval (contours), superimposed on the DSSR image (greyscale). Middle 
panels: Averaged T^in the velocity range ~ -17.5 to -15.5 km s"' (Component 2). Contour levels start at 0.7 K (23 rms) and the 
contour spacing is 1 K. Lower panels: Averaged T^in the velocity range ~ +14 to +24 km s"' (Component 3). Contour levels start 
at 0.35 K (23 rms) and the contour spacing is 0.45 K. In all cases, the molecular emission grayscale goes from 0.35 K to 3.9 K. The 
orientation of the equatorial system is indicated in the top left panel. 



A direct com parison of Component 1 with the molecular 
cloud detected by Yamaguchi et alj d 19991) (see Fig. 2u from that 
work), shows that the angular size of Component 1 is about a fac- 
tor 3-4 greater than the latter. Clearly, only the densest part of 
the molecular cloud associated with NGC 3503 (clump A) wa s 
detected in the '^CO observations of lYamaguchi et al ] (ll999h . 
From here onwards, the analysis of the molecular gas associated 
with NGC 3503 will focus on Component 1 . 



3.1.2. Kinematics and excitation conditions 



The kinematics of Component 1 was studied by using position- 
velocity maps across selected strips. The map obtained along 
b = +0°12 (corresponding to clump A) is shown in the upper 
panel of Fig. H] A noticeable velocity gradient is observed at 
velocities from about -23 km s"' to about -28 km s"' . 
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Fig. 4. Upper panel: Velocity-Galactic Lotigitude tnap obtained 
in a strip along b - +(f. 1 1 (corresponding to clump A) showing 
r^. Contour levels start at 0.7 K and the contour spacing is 0.7 
K. The dotted lines indicate the location of NGC 3503 . Lower 
panel: Mean velocity map of Component 1. Contours levels go 
from -26 to -24 km s ' , with interval of 0.2 km s"' . The dotted 
ellipse depicts approximately the region of clump A. 



Additional support to the existence of this gradient can be 
obtained through a moment analysis. Due to the large angular di- 
mensions of Component 1 and the spatial sampling of NANTEN 
observations, 28 independient CO profiles were observed to- 
wards the region enclosed by the 0.42 K-contour line in Fig. |3] 
Having these spectra a high signal-to-noise ratio (S/N > 10), 
they are well suited to study in some detail possible variations 
of the CO profiles across Component 1 . With this aim, we used 
the AIPS package to calculate the first three moments (integrated 
area, temperature weighted mean radial velocity, and velocity 
dispersion) of the CO profiles. In the lower panel of Fig. |4] the 
temperature weighted mean radial velocity distribution of CO is 
shown. A clear velocity gradient is observed along Component 
1 , and particulary in the region of clump A (depicted by the dot- 
ted ellipse), with mean radial velocities more negative toward the 
position of NGC 3503 . A velocity shift of about ~ 1.4 km s"' 
is observed across the region of clump A, which at a distance of 
2.9 kpc translates to a gradient w w 0. 15 km s"' pc"' . This panel 
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Fig. 5. Spatial distribution of Texc for Component 1 . Contours 
leves are 5.7, 7.9, 10.2, 12.4, 14.5, and 16.6 K. 



also shows that the CO emission corresponding to clump B does 
not show a significant velocity gradient. 

In order to off'er a complete picture of the kinematical prop- 
erties of clump A, we show in Fig. |6] the spatial distribution of 
the CO emission within the velocity range from -27.4 to -23.4 
km s"' . Every image depicts mean T^-values (in contours) over 
a velocity interval of 1 km s ' superimposed on the 8.13 jum 
MSX emission (in greyscale). In the velocity range from -27.4 
to -26.4 km s"' the molecular emission arising from clump A 
is slightly displaced from the brightest MSX emission located at 
(l,b) Si (289?5, +0?11). As we move towards more positive ve- 
locities, the maximum of the CO emission gradually shifts west- 
wards. 

Following lUrquhart et alJ (l2009l) . we analyze the excitation 
temperature (Texc) obtained from the CO data to probe the sur- 
face conditions of Component 1 . The excitation temperature of 
the '^CO line can be obtained considering that '^CO is optically 
thick (Ty » 1). Then, the peak temperature of this line is given 
by 



Tpeak (' CO) - Jv(Texc) - Jv(Tbg) 



(2) 



( lDickmanlll9"78l) where Jy is the Planck function at a frequency 
V. Assuming gaussian profiles for the '^CO line, combining the 
order zero moment map (i.e., integrated area) with the order 
two moment map (i.e., velocity dispersion), and using Eq. |2] 
we obtain the Tgxc distribution map (Fig.|5]l. As expected, the 
Texc distribution is quite similar to the CO emission distribu- 
tion of Component 1, reaching a maximum of ~ 17.7 K at 
(l,b) (289?45, +0?11). It is worth noting that the obtained value 
of Te^c t owards the center of c lump A is lower than the ob- 
tained bv lUrquhart etalJ (l2009h towards (l,b) ^ (289?5, -hO?11) 
(Texc = 23.9 K). This difference may be explained in terms of 
a beam smearing of our NANTEN data, which implies that the 
values of Texc shown in Fig. |5] must be considered as lower lim- 
its. 

3.1.3. Masses and densities 

The mass of the molecular gas associated with NGC 3503 can 
be derived making use of the empirical relationship between the 
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Fig. 6. Overlay of mean T^-values (contours) in the velocity range from -27.4 to -23.4 km s ' and the MSX band A emission 
(greyscale). Every image represents the CO emission distribution averaged in a velocity interval of 1 km s The velocity interval 
is indicated in the bottom right corner of each image. The lowest temperature contour is 0.84 K (~ 12 rms). The contour spacing 
temperature is 1 .4 K. 



molecular hydrogen column density, A^(H2), and the integrated 
molecular emission, /i2co 

(= / r*t/v). The conversion between 
InQQ and A^(H2) is given by the equation 



N(U2) = (1.9 + 0.3) X 10^" /nco (cm^O 



(3) 



(iDi^el et al.lll996t IStrong & Mattoxlll996h . The total molecular 
mass Mtot, was calculated through 



Mt, 



(OTj 



(Mo) (4) 



where m„,„ is the solar mass (~ 2 x lO-'"' g), // is the mean molec- 
ular weight, assumed to be equal to 2.8 a fter allowance of a rela- 
tive helium abundance of 25% by mass (lYamaguchi et al 
lYamamoto et al.ll2006l) . mn is the hydrogen atom mass (~ 1.67 
X lO"-^'* g), n is the soUd angle subtended by the CO feature 
in ster, d is the distance, expressed in cm, and Mtot is given in 
solar masses. For Component 1, we obtain a mean column den- 
sity of A^(H2) = (1.9 + 0.3) X 10^^ cm~2, and a total molecular 
mass of Mtot - (7.6 ± 2.1) x 10^ Mq . The uncertainty in M,o, (~ 
28 %) stems from the distance uncertainty and from the error 
quoted for the coefficient in Eq. [3] Areas having T"^ > 0.42 K 
were taken into account. The difference with the value derived 



by 'Yamaguch i et al.l ( Il999h . (500 Mq ), is due to both the sub- 
sampled observations of these authors and the fact that '^CO 
line is a better tracer of high density regions. 

The mean volume density (nnO of Component 1 can be de- 
rived from the ratio of its molecular mass and its volume con- 
sidering that the volume of Component 1 is the result of the 
addition of an ellipsoid with mayor and minor axis of 7'5 and 
4:5, respectively, centered at (l,b) * (289?42, H-0?15) (which in- 
cludes clump A), and a sphere of 4'5 in radius, centered at 
(l,b) ^ (289?33, +0?04) (which includes clump B). We derived 
njjo ~ 400 cm"^. Taking into account the mass and volume un- 
certainties, we derive a conservative density error of about 60 % 
(~ 240 cm-3). 

3.2. Radio continuum emission 

Figure |7] shows the radio continuum images at 843, 4800, and 
8640 MHz (in contours) superimposed on the DSSR image (in 
grayscale) in a region of ~ 9' x 9' centred on NGC 3503 . 

The continuum images show an extended source coinci- 
dent with NGC 3503 . The maxima at the three frequencies co- 
incide with the brightest optical emission region. The source, 
which has a good morphological correspondence with the op- 
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Fig. 7. Upper left panel: Radio continuum image at 843 MHz (contours) superimposed to the DSSR image (grayscale). Contours lev- 
els go from4 mJy beam"'(~ 3 rms) to 12 mJy beam^'in steps of 4 mJy beam"', and from 20 mJy beam^'in steps of 10 mJy beam"'. 
Upper right panel: Radio continuum image at 4800 MHz (contours) superimposed to the DSSR image (grayscale). Contours levels 
go from 2.4 mJy beam"'(~ 3 rms) to 10.4 mJy beam" 'in steps of 2 mJy beam"', and from 10.4 mJy beam"'in steps of 4 mJy beam"'. 
The symmetry axis of the nebula is depicted by the dotted line. Lower left panel: Radio continuum image at 8640 MHz (contours) 
superimposed to the DSSR image (grayscale). Contours levels go from 15 mJy beam"'(~ 3 rms) in steps of 1 mJy beam"'. Lower 
right panel: Emission measure profile obtained from the 4800 MHz radio continuum image along the symmetry axis of NGC 3503 . 



tical emission, exhibits a cometary shape with its symmetry 
axis indicated by the dotted line in the image at 4800 MHz 
(Fig. |7] top right panel). The fainter emission area detected 
at 843 and 4800 MHz to the northeast of NGC 3503 coincides 
with MBO and appears to be linked to the emission of the 
nebula at {l,b) ~ (289?53, +0^13). The emission region detected 
at (l,b) ^ (289?48, +0;08) at 843 MHz depicts some morpho- 
logical correspondence with the bright rim of SFO 62. The 
point-like source placed at {l,b) » (289';49,-i-0?15) without op- 
tical emission counterpart is labelled G289.49+0.15 in Fig. |2l 
From its flux densities Fg43 - 100 mJy, f48oo = 17.1 mJy, and 
^^8640 = 11.5 mJy, we derived a spectral index a = -0.92 + 
0.03 (Svocv"'), which suggests that G289.49+0.15 is a non- 
thermal extragalactic source. We substracted the emission of this 
source in order to calculate radio continuum flux densities of 
NGC 3503 . The results are given in Table [T] The spectral in- 
dex based on these estimates is a = -0. 12 + 0.03, typical for the 
thermal free-free radio continuum emission of an Hii region. 



Table 1. Radio continuum flux densities measurements of 
NGC 3503 



Frequency 


843 MHz 


4800 


MHz 


8640 MHz 


Flux density (mJy) 


320 ± 30 


270 


±50 


240 ± 40 



We estimate the emission measure EM J n^ dPj along the 
symmetry axis of the nebula using the image at 4800 MHz. The 
brightness temperature (T/,) is related to the optical depth (t) by 

ri(y) = r,x(i-o (5) 

where = 8100 + 700 K dOuireza et aT]|2006l) is tiie electron 
temperature and the optical depth t is given by 

T = 8.235 X 10"- r;'-^^ V"- ' EM (6) 
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Fig. 8. Left panel Composite image of NGC3503and its environs. Red and green show emission at 8.13 //m (MSX) and 24 fim 
(MIPSGAL). The color scale goes from 1x10"^ to 4x10"^ W m"^ sr"' and from 27 to 100 MJy ster"', respectively. The white 
rectangle encloses the IR counterpart of NGC 3503 (IRK). Right panel: Composite image of NGC 3503 . Red, green and blue show 
emission at 8 ^im (IRAC-GLIMPSE), 24 and 70 yum (MIPSGAL), respectively. Colors scales range from 40 to 500 MJy ster"' 
(24 iim), from 300 to 2000 MJy ster"' (70 //m), and from 30 to 400 MJy ster"' (8 yum). The location of the members of Pis 17 is 
indicated with black crosses. MSX and 2MASS candidate YSOs are indicated with red and green crosses (see text). 



In this last expresion, v is given in GHz, and EM in pc cm"^. 
The EM profile is shown in the lower right panel of Fig. Q As 
expected, two maxima are observed, one at the brightest sec- 
tion of NGC 3503 and the other coincident with MBO, indicat- 
ing that the electron density is higher towards these regions than 
in the sorroundings. The EM profile of NGC 3503 is consistent 
with the cometary shape morphology of the nebula seen in the 
radio continuum and optical images. The EM profile shows a 
sharp border towards the west of the nebula, while towards the 
east it decreases smoothly, indicating the existence of an elec- 
tron density gradie nt. The observed EM profile is in agreement 
with the results of ICopetti et al.l (l2000l) . who reported the ex- 
istence of an electron density gradient in the east-west direc- 
tion. Using the peak values EM - 18000 pc cm"^ (obtained at 
(l,b) « (289?50, +ail) for NGC 3503) and EM = 7000cm-'' 
(obtained at (l,b) « (289?53, +ai8) for MBO), and considering 
a pure hydrogen plasma with a depth along the line of sight 
equal to the size in the plane of the sky, electron density esti- 
mates «f = 75 + 14 cm"^ and = 45 ± 9 cm"-' are obtained for 
NGC 3503 and MBO, respectively. 

As a different approach to NGC 3503 , the rms electron 
density and i onized mass can be obtaine d using the spheri- 
cal model of Me zger & Henderson! ( Il967h and the flux den- 
sity at 4800 MHz. Assuming a constant electron density, and 
a radious /?hii = 1.25 pc, we obtain = 54 + 13 cm"-' and 
Mhii = 1 1 + 3 Mq . Thus, the electron density of NGC 3503 is 
about a factor of 5 lower than the density of Component 1 , which 
clearly indicates that the nebula has expanded. An estimate of 
the filling factor (/ = -y/«77n^) can be obtained by taking into ac- 



count the maximu m electron density d erived from optical Unes 
(«; = 154+52 cm-^: ICopetti"etani2000h . We obtain / = 0.5 - 0.8. 
Then, the ionized mass for / = 0.5 - 0.8 is in the range 8-10 
Mo. 

Regarding MBO, rms electron densities and masses es- 
timated from the image at 843 MHz are «e ^ 33 cm"-' and 
Mhii - 6 Mq , respectively. 



3.3. Infrared Emission 

The emission distribution at 8.13 pm (MSX- A band) superposed 
to the image at 24 pm (MIPSGAL) is shown in the left panel of 
Fig-E The MSX Band A includes the strong emission features at 
7.7 and 8.6 pm attributed to PAH molecules, which are consid- 
ered tracers of UV-irradiated Photodissociated Regions (PDR) 
(iHollenbach & TielensllT997! ). This image displays a small and 
strong feature, from hereonwards dubbed the IR Knot, or IRK 
for short (indicated with a white rectangle), which is coinci- 
dent with the location of NGC 3503 , and a weaker and more 
extended emission region detected in the north-western area of 
the image. The last feature will be refeiTed to as the Extended 
IR Emission, or EIE for short. IR emission at 8.3 and 24 pm ap- 
pears mixed along the whole feature. A region of low IR emis- 
sion is seen between the IRK and the EIE. The bright rim of 
SFO 62 is detected in the IR as the bright filament observed 
from {I, b) ^ (289?42, -hO?03) to (/, b) ^ (289?50, +0?08), which 
appears to be the southernmost boundary of EIE. The presence 
of PAH emission coincident with this bright rim suggests the ex- 
istence of a PDR at the southern edge of Component 1 . 
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Fig. 9. Upper panel: distribution estimated from HIRES 
IRAS images at 60 and 100 jum. Contour levels go from 25 to 
31 K in steps of 1 K, and from 33 K in steps of 2 K. Lower 
panel: Overlay of the distribution (contours) and the DSSR 
image. 

The right panel of Fig. [8] shows a detailed image of the IRK, 
at 8.0 yum {Spitzer-IRAC band 4, in red), 24 yum (MIPSGAL, 
in green), and 70 jdva (MIPSGAL, in blue). The emissions at 8 
and 70 //m show a bright half shell-like feature, which encloses 
the position of the stars of Pis 17, indicated as black crosses. 
Arc-shaped faint filaments are also detected in the 8 jim emis- 
sion in the northeastern section of IRK. The morphology of 
the IR emission is compatible with the cometary-shape of the 
Hii region and very likely indicates the presence of a PDR be- 
tween NGC 3503 and clump A. The MIPSGAL emission at 24 
yum, which is projected at the center of the structure, indicates 
the existence of warm dust inside the Hii region (as the cases of 
NIO and N2 1 . IWatson etal]|2008l) . 

IRAS 60 and 100 pm data show dust with color tempera- 
tures between about 20 to 100 K, which correspo nds to the "cool 
dust component" (see lSreenilavam & Fichll201 1 , and references 
therein). A dust temperature (Td) map was produced using the 
equation 



Td = 95.94 / ln(B„) (K) 



(7) 



dPraind 119901: IWhittetl 119921; ICichowolski et al.l 1200 lb . where 
B„ - 1.667^^^"'(Fi()o/F6o) is the modified Planck function, with 
Fioo and Feo being the 100 pm and 60 pva fluxes, respectively. 
The parameter n = 1.5 is related to the absorption efficiency 
of the dust {ky oc v", normalized to 40 cm^ g"' at 100 pm). 



The obtained dust temperature map is shown in the upper panel 
of Fig. |9] Dust temperature goes from 25 K to 46 K. Fig. |9] 
(lower panel) also shows that the region with the highest dust 
temperatures (~ 46 K) coincides with the brightest section of 
NGC 3503 . This temperature is in go od agreement wit h those 
obtained in RCW 121 and RC W 122 jArnal et alj|2008h and in 
NGC 6357 (Capp a et al]l201lh although they are slightly higher 
than typical valu es for Hn r egions ( ~ 30 K) (see Cappa et al] 
2008; Cichowolski et al.l|2009; Vasgu ez et alj|2010h . The stellar 
UV radiation field of the stars in NGC 3503 is responsible for 
the heating of the dust immersed in the Hii region. A noticeable 
feature is a fringe of cold (~ 25 K) dust suiTounding the po- 
sition of NGC 3503 , which coincides with the region showing 
low emission at 8.3 pm placed between IRK and EIE described 
above. An extended region of warmer (~ 36 K) dust is detected 
to the northeast of NGC 3503 , coincident with the position of 
MBO, indicating that this optical feature has a radiatively heated 
dust component. In Table |2] we summarized the IR and dust pa- 
rameters for IRK and EIE. The averaged dust temperatures (Td) 
are calculated using Eq. Qand flux densities at 60 and 100 pm 
obtained integrating several polygons around the sources. Dust 
masses are derived from 



Md = m„ Ffio d' (Bf - 1) (Mq) 



(8) 



where d the distance in kpc, Feo is given Jy, and 
OTi.5 = 0.3 X 10"^. Based on the molecular gas mass derived for 
Component 1, the ionized gas of NGC 3503 and MBO, and the 
dust masses obtained for the IRK and the EIE (~ 10 M© ), we 
derive a mean weighted gas-to-dust ratio of about ~ 800 + 300. 
This ratio is a factor of about ~ 5 higher than the average values 
assumed for the Galaxy (140 + 50; iTokunaga.2000.) . Note how- 
ever that dust with temperatures < 20 K ("cold dust compo- 
nent"), which is detected at wavelengths not surveyed by IRAS 
{A> 100 pm), dominates the dust mass by a factor larger than ~ 
70 with respect to the cool dust component (ISreenilavam & FichI 
i201 1,) . This notoriously decreases the gas-to-dust ratio. The dust 
temperature estimated for EIE may suggest the existence of ion- 
izing sources inside Component 1. However, a search for OB 
stars embedded in Component 1 performed using the available 
VIZIER catalogues failed to detect such sources. 

To investigate the presence of protostellar candidates possi- 
bly related to NGC 3503, we used data from the MSX, 2 MASS, 
and IRAS point source catalogs. We searched for point sources 
in a region of about 10' in size centered at the position of the Hii 
region. Taking into account sources with flux quality ^ > 2, a to- 
tal of 69 MSX point sources were found projected onto the area. 
Based on F2i/Fg, and F14/F12 ratios, where Fg, F12, Fu, and F21 
are the fluxes at 8.3, 12, 14, and 21 ^ m, respectively, and af - 
ter applying the criteria summarized bv lLumsden et al.l (l2002h . 
we were left with three sources with F2\IF% > 2 and f 14/^12 < 
1, which are classified as compact Hn regions (CHn). This is 
indicative of stellar formation. The location of these sources is 
indicated with red crosses in Fig.jSl 

Using the 2 MASS catalog (iCutri et al.1 l2003h . which pro- 
vides detections in J, H and Ks bands, we searched for 
point sources with infrared excess. Taking into account sources 
with signal-to-noise ratio (S/N) > 10 (corresponding to quality 
"AAA"), we found 5528 s ources projected onto a circular region 
of 5' in radius. Following lComeron et al ] d2005h . we determined 
the parameter q- { J - H) - 1.83 x (i/ - Ks). Sources with q < 
-0.15 are classified as objects with infrared excess , i.e. can- 
didate YSOs. By applying the above criteria, we found only 9 
sources projected onto NGC 3503 . The location of these sources 
is also indicated in Fig. |8] with green crosses. The location of the 
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Table 3. Candidate YSOs obtained from the MSX and 2 MASS catalogues. 





(hb) 


MSX source 


Fs (Jy) 


Fn (Jy) 


F,4 (Jy) 


F21 (Jy) 


1 


289.48,+0.14 


G289.4859+00.1420 


1.49250 


1.8145 


0.9005 


0.3180 


2 


289.49,+0.12 


G289.4993+00.1231 


0.82699 


1.4887 


1.2237 


4.3818 


3 


289.50,+0.11 


G289.505 1+00. 1161 


0.87911 


2.8632 


1.6759 


5.6330 




iUb) 


2 MASS source 


J (mag) 


H (mag) 


K (mag) 




4 


289.49,+0.14 


11011934-5949422 


13.277 


13.047 


12.806 




5 


289.50,+0.15 


11012021-5949123 


15.084 


14.837 


14.506 




6 


289.50,+0.13 


11011933-5950037 


14.539 


14.100 


13.696 




7 


289.50,+0.14 


11012329-5949269 


14.862 


14.380 


14.027 




8 


289.54,+0.13 


11013591-5951037 


15.237 


14.592 


14.106 




9 


289.53,+0.12 


11013263-5951140 


15.633 


14.946 


14.404 




10 


289.52,+0.10 


11012238-5952449 


12.747 


11.744 


10.836 




11 


289.52,+0.15 


11013238-5949339 


14.754 


14.243 


13.873 




12 


289.48,4-0.14 


11011390-5949099 


14.888 


13.661 


12.827 





Table 2. Main infrared parameters inferred from IRAS fluxes at 
60 yum and 100 /im 



Parameter 



IRK 



EIE 



Sso (Jy) 
Sioo (Jy) 

frf(K) 
Mrf(Me) 



1470 ± 20 700 ± 15 

2170 ± 25 1970 ± 25 

-37 -28 

-3 -7 



2 MASS source 11011390-5949099 (source #12) is almost co- 
incident with the MSX source G289.4859-H00. 1420 (source #1). 
These sources are projected onto an intense mid-IR clump lo- 
cated in the northeastern border of the bright half shell-like fea- 
ture, at (/, b) = (289?48,-H0n4 ) (see Fig.|8] right panel), and onto 
a small enhancement in the radio continuum emission detected 
at 4800 MHz (see Fig. [TOl below). These characteristics make 
this object an excellent candidate for investigating star forma- 
tion with high angular resolution observations. The names of the 
YSO candidates, their position, fluxes and magnitudes at differ- 
ent IR wavelengths are listed in Table [3] 

4. Discussion 

4.1. The ionized gas 

The number of ionizing Lyman continuum photons (A^Lyc) 
needed to sustain the ioniz ation in NGC 3503 can be calculated 
using the equation given bv lSimpson & Rubinl(ll990l) 



A^Lyc = 4.45 X 10-*** 5(4800 MHz) 



(9) 



(ISternberg et alJl2003h . which is capable of ionizing NGC 3503 
in a greement with previ ous results by [Thompson et al.l (l2004t) 
and iPinheiro et alj ( 1201 0). The last authors asserted one BO V 



star and three B2 V stars, belonging to the open cluster Pisl7, 
inside the nebula as ionizing sources. Since the ionizing photons 
emitted by a B2 V star are significantly fewer than those of a BO 
V star, their contribution to the energetics of the nebula can be 
neglected. 

In Sec. 3.2 we pointed out on the discrepancy between the 
electron and molecular densities. The density of Component 1 is 
about a factor of 5 higher than the electron density of NGC 3503 , 
which clearly indicates that NGC 3503 has been expanding as 
a result of unbalanced pressure between ionised and molecular 
gas. In order to estimat e the dynamical age (fdvn) of the Hii region 
we used the model of iDvson & Williamsl (Il997h . The radious of 
an Hii region (Rmi) in a uniform medium is given by 



4/7 



(10) 



were Rs is the radious of the Stromgren sphere (IStromgrenll939h 

/ , \l/3 

before expanding, given by R^ - \liN*^^,^IAn {2n^^Yap\ , and 

Vs is the sound speed in the ionized gas (~ lOkms '). For 
Rmi - 1-25 pc (1.'5 at a distance of 2.9 kpc), an ambient den- 
sity «H, = 400 + 240 cm"^ (see Section 3.1.3), and A^^^.^ = ^-^1 
X 10"*^ s we infer fdyn - 2 x 10^ yr. The expansion velocity of 
the Hii region can be estimated by means of 7?hii = v , ( I , 

yielding an expansion velocity of about - 5 km s The ob- 
tained dynamical age and expansion velocity are i n agreement 
with those obta ined in typical Hii region s (Gum 31. ICappa et al.l 
[2008: Sh2-173. lcichowolski et al.ll2009h . 

The total number of ionizing Lyman continuum 
photons needed to sustain the ionization in MBO is 



Lyc 



where d is the distance in kpc, and 5(4800 mhz) is the flux den- 
sity at 4800 MHz in Jy. We obtain A^Lyc = (1.8 + 0.4) x 10^^ s"' . 
This number is a lower limit to the total number of Lyman 
continuum photons required to maintain the gas ionized, since 
about 25 - 50 % of the UV photons are absorbed by inter- 
stellar dust in the Hii region (Inoue 2001). Consequently, we 

need A^Lyc - 3.6 x 10 s"' . The number of Lyman continuum 

photons emitted by a BO V star is A^*^^ = 1.07 x lO"*** s^' the star are A?*J^ = 3.2 x 10'*^ s^' (ISmith et al.ll2002l) . Although 



2.5 X 10 seg . To search for stars that can pro- 
vide the necessary UV photons, we used the available 
VIZIER catalogues. CP-5 9 2951 is a Bl V star p laced at 
f/,fe) = (289°511, H-0?195) (iBuscombe & Fosteii [T999I) . Using 
the catalogued spectral type, the v isual magnitude, and the 
calibration of Schmidt-Kale^ (Il982h . we estimated a distance 
d = 2.7 ± 0.8 kpc. The continuum ionizing photons emitted by 



Lyc 
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the distance of this star is in agreement (within errors) with the 
distance of NGC 3503 , the value of A^^^^ is almost one order of 
magnitude lower than the required to ionize MBO. The presence 
of Pisl7 at a projected distance of 3.4 pc (~ 4' at a distance of 
2.9 kpc) indicates that the contribution of the star cluster to the 
ionization of MBO can not be ruled out. 



4.2. Star formation process 

The presence of 12 candidate YSOs projected onto the 
NGC 3503 -region suggests that they may have been triggered 
by the expansion of the Hii region through the "collect and col- 
lapse" model, which indicates that expanding nebulae compress 
gas between the ionization and the shock fronts, leading to the 
formation of molecular cores whe re new stars can be embed- 
ded. Using the analytical model of IWhitworth et all (11994) for 
the case of expanding Hii regions, we derived the time when 
the fragmentation may have occurred (t/rag), and the size of the 
Hn region at tf,-ag (Rfrag), which are given by 



Rfrag[pc] 



(11) 



(12) 



where 02 is the sound velocity in units of 0.2 km s , = 
uhJIOOQ, and N49 = Nl^JlQ*"^. Adopting for this region 0.3 km 

s~^ for the sound velocity, which corresponds to temperatures 
of 10-15 K in the surrounding molecular clouds (see Section 
3.1.2), we obtained f/ra^ ~ 3.5 x 10^ yr, and Rfrag ~'!.5pc. 
Considering that the values of tfrag and Rfrag are larger than tdvn 
and Rhh (see Section 4.1), we can conclude that the fragmenta- 
tion in the edge of NGC 3503 is doubtful. 

4.3. Kinematics of clump A 

In section 3.2 we reported the existence of a velocity gradient 
across clump A. Velocity gradients in molecular cores/clumps 
were usually interpreted as gravitationally bound rotation mo- 
tions (particulary in dense and small molecular cores). Several 
theoretical models predict cloud flattening perpendicular to the 
rotation axis in response to centrifugal stress, which at first 
glance seems to be suitable for clump A. In order to investi- 
gate the dyn amical stability of clum p A, we use the parameter 
jS defined by iGoodman et al.l (Il993l) to quantify the dynamical 
role of rotation by comparing the rotational kinetic energy to the 
gravitational energy. Thus, /3 can be written as 



(1/2) I oj'^ 
qGM^/R 



(1/2) (p/q) 



R' 



GM 



(13) 



where / is the moment of inertia (/ - pMR^), qOM^/R is 
the gravitational potential energy, and co' - a>/sin(i), where 
i is the inclination of the cloud along the lin e of sight. 
Considering p/q-0.22 (see iGoodman et al.l [T993h . sin(0 = 1, 
w=0.15 km s~' pc'' (see Sect 3.1.2), /J ~ 4.5 pc, and a lower 
Hmit mass to clump A (M ~ 3 x 10^ Mq ), we obtain /3 « 0.02. 
This extremely low value of /? indicates that the effect of rotation, 
if exists, is not significant in mantaining the dynamical stability 
of clump A. This might weaken the rotating cloud interpreta- 
tion. Furthermore, Figs.|4]and|6]show that only molecular gas at 

more negative velocities ( 27 km s"' / -26 km s"' ) coincide 

with NGC 3503 , which suggests that the velocity gradient might 
be a direct consequence of an interaction between clump A and 
NGC 3503. 



Therefore, although rotation can not be entirely ruled out 
with the present data, we consider different origins for the ve- 
locity gradient observed across clump A, namely, 1) the expan- 
sion of the nebula at ~ 5 km s"' (see Sect 4.1) has been ac- 
cumulating molecular gas behind the shock front which orig- 
inated the expansion of clump A at approximately the same 
velocity of the ionized gas, as expected according to the mod- 
els of Hosokawa & Inutsuka (2006), 2) the velocity gradient of 
clump A is the consequence of a collision between Component 
1 and another molecular cloud, which in turn, might have in- 
duced the formation of Pis 17 and the candidate YSOs reported 
in Sect. 3.3. Although they are rare, cloud-cloud collisions can 
lead to gravitational instabilities in th e dense, shocked gas, re- 
sulting in triggered star formation (see lElmegreenI 1 998l and ref- 
erences therein), 3) clump A is actually composed by different 
subclumps at different velocities, which are not resolved by the 
NANTEN observations. In this case, NGC 3503 might be related 
with a subclump at more negative velocities. 

Further high-resolution studies with instruments hke APEX 
may help to clarify this question. 

4.4. Tlie PDR at tlie interface between NGC 3503 and clump 
A 

As mentioned in Sect 3.1.2, clump A exhibits an excita- 
tion temperature Tg^-c ^ 17.7 K, which is also the highest ex- 
citation temperature along Component 1. This temperature 
is higher than expected inside molecular cores if only cos- 
mic ray ioniz ation is considered as the m ain heating source 
(r ~ 8 - 10 K. lvan der Tak & van Disho eck 2000). which im- 
plies that additional heating processes are present close to clump 
A. Very likely, clump A is being exte rnally heated through the 
photoionisation of its surface layers (lUrquhart et al.l l2009l) as 
a consequence of its proximity to NGC 3503 . This scenario is 
in line with the presence of the PDR at the interface between 
NGC 3503 and clump A (see Sect. 3.3). In this context, it would 
be instructive to make a simple comparison of the PDR dust sur- 
face temperature, with the dust temperature obtained before to- 
wards the edge of NGC 3503 (see Fig©. 

The structure of a PDR is governed by the intensity of the 
UV radiation field impinging onto the cloud surface (Go) in units 
of the Habins ( 1963) F UV flux (1.6 xlO'^ ergs cm"^ s '), with 
Go <x N^^l^ nT y (lTielensl2005h . where is the fraction lu- 
minosity over 6 eV, and is the luminosity of the star Ad opting 
A^Ly = 1.07x10^^ s-1 and L = 7.6 x 10"* Lo (Sternberietay 
l2003h . X - and considering «g ^ 75 - 154 cm""* (see Sect. 
3.2), a range of Go - (0.5 - 1.5) x 10"^ is obtained. The dis- 
tribution of dust temperature in the PDR (T^'^'^) is governed by 
the absorption of the stellar photons which are reemited by dust 
grains as IR photons. Following Tielens (.2005 ). we can relate 
r''"'' with the visual absorption along the PDR (Ay) as 



(T'fy - vo Go e^-'-''^^' + ln(3.5 x lO-^rjoo To) Turn To" (14) 

where vo is the frequency at 0.1 fim, tioo = lO'-' is the effec- 
tive optical depth at 100 //m, and To is the temperature of the 
slab estimated as To = 12.2 Go''^. Considering both, the posi- 
tion of NGC 3503 at the edge of clump A, and that the low value 
of visual absorption derived for the members of Pis 17 (Ay = 
1.6 mag. iPinheiro et al] 12010 *) is probably due to the interstel- 
lar absorption in the line of sight of the Hii region, we can as- 
sume Av ~ for the surface of the PDR. Taking into account 
the range of Go and Eq. [141 we obtain T'^'^' ^40-55 K, in ac- 
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cordance with the observed dust color temperature estimated for 
NGC 3503 (~46 K). This gives additional support to the exis- 
tence of a PDR between NGC 3503 and clump A. 

4.5. Possible scenario 

Figure [10] displays a composite image of NGC3503 and its en- 
virons. White and light blue contours show the radio continuum 
emission at 4800 MHz coiTesponding to NGC3503 and MBO, 
and the CO emission, respectively. The color scale shows the 
emission at 8 //m (red) and 4.5 fim (green). 

As described in Sect. 3.3 and 4.3, emission attributed 
to PAHs encircles the southern and western borders of the 
Hn region, indicating the location of the PDR. The position 
of NGC 3503 near the border of clump A and the existence 
of an electron density gradient along the symmetry axis of 
the Hi: region (with the higher electron densities close to the 
strongest CO emission) suggest that the ionised gas is being 
streamed away from the densest part of the molecular clou d. This 
is in dicative that NGC 3503 is a blister-type Hn region (llsraell 
11978) that probably has undergone a champagne phase. 

The so-called Champagne flow m odel (Tenorio-Tagl3 197^ 
iBediin & Tenorio-Tadei [l98lt iTenorio-Tagle & BediinI (l982l) 
proposes that an expanding Hii region placed at the edge of a 
molecular cloud eventually reaches the border of the cloud and 
expands freely in the lower density surrounding gas, originat- 
ing an extended Hn region with a characteristic density distri- 
bution. The Hn region becomes density bounded towards the 
lower density region, while it is ionization bounded towards 
the molecula r clou d. This scenario was formerly proposed by 
ICopetti et al.l (12000) to explain the electron density gradient ob- 
served across NGC 3503. The location of Pis 17 close to the 
brightest radio continuum region is consistent with a projec- 
tion effect ( Yorke et al. 1983). In this scenario. Pis 17 originated 
NGC 3503, which reached the northeastern border of clump A 
after 2x 10^ yr This time represents a lower limit to the age 
of the Hn region, since the inferred main-seque nce life time o f 
the main ionizing star is about (3 -5) x 10^ yr ('M assevl[l998h . 
The leakage of ionized gas and UV photons might have con- 
tributed to the formation of MBO, since its location (along the 
symmetry axis of NGC 3503 ) and its low density (about half of 
NGC 3503 ) suggest that this feature may consist of ionized gas 
which has scaped from NGC 3503 after a time of 2x 10^ yr. 

The velocities of both the molecular and ionized gas are 
compatible with the champagne scenario. The mean velocity of 
Component 1 (-24.7 km s"' , see Sect. 3.1.1) corresponds to the 
natal cloud where NGC 3503 originated. Molecular gas having 
more negative velocities is mainly linked to clump A and may 
represent material moving away from the ionized gas placed in 
front of the Hn region (as a result of expanding motions or an 
external impact over Componen t 1). Th e ionized gas, having a 
velocity of -21 km s"' (Georgel in et al.l l2b00). might be reced- 
ing the ob server. We note, however, that the Fabry-Perot obser- 
vations bv lGeorgelin et al ] (12000 ) have a spectral resolution of 5 
km s"' . 

NGC 3503 resembles the two well studied cases of S305 and 
S307. These Hn regions show a non spherical morphology in the 
1465 MHz map of Fichl(ll993h . and they are placed close to '^CO 
emission peaks, which suggests that they are locate d close to the 
hottest part of their parental molecular clouds (see lRusseil et alj 
11995 . and references therein). Both objects show significant 
electron density dependence on position, which can also be in- 
terpreted as radial gradients (ICopetti et al. 2000). Furthermore, 
both Hn regions show discrepancies of about ~ 12kms ' be- 



tween the velocities of the ionized and molecular gas, prob- 
ably due to a champagne effect producing a flow toward the 
observer (Russeil et al. 1995). In the case of S307, the cham- 
pagne scenario is reinforced by its half-she ll shape, and the high 
density in the brig htest part of the nebula dFelli & Harten|[T98ll : 
lAlbert et alJ[T98a) 

According to the champagne scenario, a velocity gradi- 
ent along the sym metry axis of the nebula can be expected 
dGarav et alj|1994l) . In the case of NGC 3503 increasing veloci- 
ties with the distance to the molecular cloud are expected. Thus, 
high spectral resolution radio recombination are required to an- 
alyze the ionized gas velocities and may give additional support 
to our interpretation. 

5. Summary 

NGC 3503 is a bright Hii region of ~ 3' in size centered at 
(l,b) = (289';51, -h0';12) located at a distance of 2.9 + 0.4 kpc. 
The ionizing sources are B-stars belonging to the open cluster 
Pis 17. 

With the aim of investigating the molecular gas and dust 
distribution in the environs of the Hii region and analyzing 
the interaction of the nebula and Pis 17 with their molecular 
environment, we analyzed '^CO(l-O) data of a region of 0?7 
in size obtained with the NANTEN telescope with an angular 
resolution of 2.'7, radio continuum data of NGC 3503 at 4800 
and 8640 MHz obtained with the ATCA telescope (with syn- 
thesized beams of 21" and 13", respectively), and data at 843 
MHz retrieved from SUMSS, and available IRAS, MSX, IRAC- 
GLIMPSE, and MIPSGAL images. 

The analysis of the CO data allowed the molecular gas linked 
to the nebula to be mapped. This molecular gas (Component 1) 
has a mean velocity of -24.7 km s"' , a total mass and density 
of (7.6 + 2.1) X 10^^ Moand 400 + 240 cm"^ respectively, and 
displays two clumps centered at (l,b) - (289^47, -i-0?12) (clump 
A) and (l,b) = (289?32, -hO?03) (clump B). The morphological 
correspondence of the molecular emission with a large patch 
of high optical absorption adjacent to NGC 3503 is excellent. 
NGC 3503 is projected near the border of clump A, with the 
strongest molecular emission adjacent to the highest electron 
density regions. The agreement of the molecular velocities with 
the velocity of the ionizad gas (-21 km s"' ) as well as the mor- 
phological correspondence with the nebula indicate that clump 
A is associated with NGC 3503. The more negative velocities of 
the gas in clump A, coincident with the Hn region, are probably 
due to the expansion of the Hn region or an external impact over 
Component 1. 

The analysis of the radio continuum images confirms the 
electro n density gradient previously found by ICopetti et al.l 
(120001) . The images show that NGC 3503 exhibits a cometary 
morphology, with the higher density region near the maximum 
of clump A. The rms electron density and the ionized mass 
amount to 54 + 13 cm"-' and 9 ± 3 Mq , respectively. A low den- 
sity ionized region (MBO) located close to the lower electron 
density area of NGC 3503 is also identified in these images. 

Strong emission at 8 /jm surrounds the bright radio contin- 
uum region of the nebula, indicating the presence of a photodis- 
sociated region at the interface between the ionized region and 
clump A. MIPSGAL emission at 24 /im shows the existence of 
warm dust inside the Hn region. Based on high resolution IRAS 
images at 60 and 100 jum, a mean dust color temperature and dust 
mass of 37 K and 3 M© were estimated for NGC 3503 . The pres- 
ence of candidate YSOs projected onto the Hll region, detected 
using MSX and 2MASS point source catalogues, suggests the 
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Fig. 10. Composite image of NGC 3503 and its environs. Red and green show emission at 8 and 4.5 fim (IRAC-GLIMPSE), 
respectively. White and light blue contours show the radiocontinuum 4800 MHz and CO line emission, respectively. 



existence of protostellar objects in the neighbourhood of NGC 
3503, although there is no clear evidence of triggered star for- 
mation. 

The location of NGC 3503 at the edge of a molecular clump 
and the electron density gradient in the Hn region suggests 
that NGC 3503 is a blister-type Hii region that has undergone 
a champagne phase. In this scenario, the massive stars of Pis 17 
created NGC 3503, which reached the border of the molecular 
cloud after 2x10^ yr. The leakage of ionized gas and UV pho- 
tons have probably contributed in the formation of MBO. Thus, 
the spatial distribution of the molecular gas and PAHs give ad- 
ditional support to a scenario first proposed by Copetti et al] 
(1200 0). The proposed scenario for NGC 3503 may explain the 
slight difference between the main velocity of its parental molec- 
ular cloud (-24.7 km s"' ) and the velocity of its ionized gas 
(-21 km s ' ). High resolution radio recombination line obser- 
vations may help to confirm the proposed scenario. 
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